Hallmark of strong electronic correlations in LaNiOa: 
photoemission kink and broadening of fully occupied bands 
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Recent angular-resolved photoemission experiments on LaNiOs reported a renormalization of the 
Fermi velocity of Cg quasiparticles, a kink in their dispersion at —0.2 eV and a large broadening and 
weakened dispersion of the occupied t2g states. We show here that all these features result from elec- 
tronic correlations and are quantitatively reproduced by calculations combining density-functional 
theory and dynamical mean-field theory. The importance and general relevance of correlation effects 
in filled bands coupled by inter-orbital interactions to a partially-filled band are pointed out. 

PACS numbers: 71.27.+a,71.15.Mb,72.80.Ga 






I 

o 
o 



> 

o 

in 

o 



% 



The effect of strong correlations in partially-filled nar- 
row bands is a central topic in condensed matter physics. 
Experimental and theoretical studies have revealed their 
key role in metal-insulator transitions, high- Tc supercon- 
ductivity and heavy-fermion behavior. The remarkable 
development of sophisticated experimental techniques of- 
fers novel insights, validates established concepts but 
also raises new questions. Recently, Eguchi et. al [l| 
performed a delicate soft X-ray angular-resolved photoe- 
mission (ARPES) study of LaNiOa , a three-dimensional 
correlated oxide. They found clear signatures of strong 
correlations in the partially filled quasiparticle band but 
also, surprisingly, in completely filled bands. 

Among the remarkable family of charge-transfer in- 
sulators RNiOa (with R a rare-earth element) [2|, |3| , 
LaNiOs is the only one which remains metallic down to 
the lowest temperatures. Renewed interest in the nick- 
elates has been triggered by the investigation of oxide 
heterostructures ^\B,, and by LaNiOs being a possible 
electron analogue of cuprate superconductors J8l4llj . The 
nominal electronic configuration of the nickel orbitals is 
^206^. First-principle calculations based on density func- 
tional theory (DFT) show that the t2g bands are fully 
occupied and that the two nearly-degenerate Cg bands 
cross the Fermi level, resulting in a two-sheet Fermi sur- 
face [Ij, ll3[ . Clear experimental evidence for strong elec- 
tronic correlations, such as enhanced effective mass and 
reduced Drude weight, are revealed by the measurement 
of the low-temperature specific- heat, magnetic suscep- 
tibility, resistivity and thermopower jl4Ml7[ as well as 



optical spectroscopy 18|-|20[. 

The ARPES experiment of Eguchi et al. provides im- 
portant additional insight by revealing several interesting 
features. First, there is a renormalization of the eg quasi- 
particles Fermi velocity (by a factor of ~ 1/3 as compared 
to the DFT band-structure) at the Fermi-surface crossing 



near the F-point. Second, the dispersion of these quasi- 
particles near F displays a 'kink' at a binding energy of 
order —0.2 eV. Finally, in the energy range where DFT 
predicts dispersive t2g bands, the ARPES spectra dis- 
play instead a broad signal which disperses weakly and 
is shifted towards larger binding energy. While the Fermi 
velocity renormalization is consistent with the enhance- 
ments observed in thermodynamics and transport mea- 
surements, these findings raise the following questions, 
(i) What is the physical origin of this renormalization 
and, especially, of the kink structure? (ii) Is the ARPES 
signal in the [—1.0, —0.5] eV range indeed associated with 
the t2g bands, and if so, why is it broad and weakly dis- 
persing? 

In this letter, we address these questions using a combi- 
nation of density functional theory and dynamical mean- 
field theory (DFT+DMFT), as well as model studies. 
Our results yield a Fermi velocity renormalization and 
a kink structure in excellent agreement with ARPES, 
which establishes that both features are a consequence 
of electronic correlations. We also find that the fully 
filled t2g bands have a weakened dispersion and undergo 
significant broadening. This challenges the naive picture 
according to which interactions only have a small effect 
on fully occupied bands. By investigating a simplified 
model, we explain this finding by the coupling of holes in 
the fully occupied band to particle-hole excitations in the 
partially occupied one, due to inter-orbital interactions. 
This observation is relevant to other correlated materials 
in which the energy separation between partially-filled 
and fully occupied bands is small. 

Our calculations use the full potential implementa- 
tion of DFT-hDMFTdetailed in Ref. iH. Well-localized 
Wannier functions 
ergy window 



[22| are constructed from a large en- 
-7.8, -1-3.2] eV, which includes all nickel 
d-states and oxygen p-states, as appropriate for such 
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FIG. 1. (Color online) Spectral density map close along T — X 
(close-up), projected onto Sg character. Dotted line (red) 
quasiparticle dispersion setting ImE — 0"*" . Plain line (green) 
DFT-GGA band. Open circles (blue): ARPES data [I[. Inset 
Real (black, thin line) and imaginary part (red, thick line) of 
the eg-orbital self-energies. Dashed lines: linear (parabolic) 
fit of the real (imaginary) part. 



a charge-transfer compound with significant hybridiza- 
tion between those states. For simplicity, most calcu- 
lations were performed for the cubic structure, ignoring 
the small rhomboedral distortion of LaNiOs . Local in- 
teractions between Ni-Sd states are included in the form 
1/2 U^^ihj-nuhm'a'- We checked that for the problem at 
hand the omitted spin-flip and pair-hopping terms are 
not important. The reduced interaction matrices U^^, 
and U^^i for equal and opposite spins respectively are 
constructed from three Slater integrals F", F^, F'^ in 
the usual manner. We use the standard approximation 
pA, lp'2 ^ 0.625 and define the interaction parameters 
U and J through U = F°, J = (F^ + F'^yU. Stan- 
dard values for nickel compounds [1^ [2j| U = 8.0 eV 
and J = 1.0 eV are used. The double-counting correc- 
tion term is taken in the 'around mean-field' form as 
customary for metallic compounds E^^ = UN — [U + 
[M — l)J]Na/M , where M is the number of correlated 
orbitals, N is the total occupancy and N^ is total oc- 
cupancy of spin a. The calculations are fully charge 
self-consistent, which turns out to be essential in or- 
der to insure that the oxygen states are positioned cor- 
rectly in agreement with PES spectra [2J]. We solve 
the DMFT quantum impurity (embedded atom) prob- 
lem using the TRIQS |25| toolkit and its implementation 
of the hybridization expansion continuous-time quantum 
Monte-Carlo algorithm [26|, |27| , using Legendre polyno- 
mials 128[ . The imaginary-time data is continued analyt- 
ically using the stochastic maximum entropy method [29[ 
and Fade approximants [30J. All calculations were done 
at a temperature T = 1/50 eV~ 230 K. 

We now describe our results. Focusing first on the 



FIG. 2. (Golor online) Spectral intensity map (all orbitals). 
DFT-GGA bands are shown for comparison (green, plain 
lines). The broadening, shifting and weaker dispersion of the 
t2g states described in the text are clearly apparent. The 
boundaries of Fig. [T] are also shown (blue, dotted). 



Cg bands, wc display on Fig. [T] the spectral intensity in 
the vicinity of the F point. For clarity, the t2g bands 
have been projected out. The quasiparticle dispersion 
obtained by setting the imaginary part of the self-energy 
to 0+ is also shown as a guide to the eye (dotted line). 
The Fermi velocity is renormalizcd, corresponding to an 
effective mass m* /m ~ 1/ Z ~ 3, in good agreement with 
experiments. As one moves further away from the Fermi 
surface crossing, a 'kink' in the dispersion is found, at an 
energy close to —0.2 eV, followed by a "waterfall" struc- 
ture corresponding to a rapid broadening of the spectral 
function. These results agree well with the X-ray ARPES 
experiments [l[, providing strong evidence for a purely 
electronic origin of these structures. Indeed, kinks orig- 
inating from strong electron correlations have been doc- 
umented in previous DMFT studies at the model level 
and for other materials J3lM33j . The kink can be directly 



related to the eg self-energies displayed in the inset of 
Fig. [TJ We note that the low-frequency Fermi-liquid be- 
havior applies only to a narrow frequency range around 
cj = 0. A clear deviation, from linearity for ReS and from 
quadratic behavior for ImE, is observed at a negative fre- 
quency Lo ^ —0.2 eV. Below this energy scale, ReE has a 
kink and bends downwards, while ImE becomes large so 
that a quasiparticle description no longer applies. 

Now we turn to the t2g bands. A full map of the 
spectral intensity including all orbitals is displayed on 
Fig. [21 for binding energies between —1 and -1-0.5 eV. 
It is seen that, in comparison to the DFT-GGA band 
structure (solid lines), the t2g states (i) are pushed fur- 
ther below the Fermi level (with e.g. the top of the band 
at ^ —0.4 eV at the F-point in contrast to ~ —0.2 eV 
in DFT-GGA ) (ii) have a weaker dispersion than the 
band structure result and (iii) importantly, undergo a 
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FIG. 3. Second-order contribution to the self-energy. 



considerable broadening, due to a rather large scatter- 
ing rate ImE^^ (w) in the energy range from uj ^ ~1 eV 
to ~ —0.5 eV, even though the t2g states are fully oc- 
cupied. All three facts agree with the ARPES experi- 
ments jl| , which report broad and weakly dispersive spec- 
tral weight along the T-X and M-R directions, and no 
evidence for dispersive i2g bands along M-R in the en- 
ergy range [—0.5, —0.2] eV where DFT-GGA would place 
these bands. A refined quantitative positioning of the t2g 
states would require to take the rhomboedral distortion 
of LaNiOa into account. 

Hence, it is seen that the effects of interactions on one- 
particle spectra are strong for the t2g bands, in spite of 
those bands being completely filled. This may seem sur- 
prising if one follows naive intuition drawn from single- 
band models. However, in a multi-band context this is 
quite natural: a hole created in the occupied t2g band (as 
in a photoemission process) can scatter against particle- 
hole excitations associated with the partially occupied 
Cg states, due to the inter-orbital matrix elements. Con- 
sider for example the second-order diagram for the t2g 
self-energy shown in Fig.[3]in which a hole in an occupied 
band /? couples through Ua^ (q) to a particle- hole excita- 
tion of band a described by the corresponding polariza- 
tion bubble. The latter is proportional to ria^k'+q — na.k' 
and would vanish for a fully filled or empty band, but 
gives a non-zero contribution if a is partially occupied. 
This conclusion persists for higher-order diagrams. This 
effect has similarities with the physics of X-ray edge sin- 
gularities and core-hole photoemission, with the impor- 
tant difference that the hole is not static in the present 
case. 

In order to understand this important effect in 
a simpler setting, we consider a two-band Hub- 
bard model H 
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<^kcl^„Ckraa ~ A E^ n,,2a + 

In this expression, m ~ 
1,2 is a band index and A a crystal-field splitting such 
that the second band is lower in energy. We take both 
bands to have a semicircular density of states with the 
same half-bandwidth D, U/D = 2.0, and no inter-band 
hybridization. We set the chemical potential so that 
there are 3 electrons per site. This model captures 
the transition from a two-band metal into a single-band 
metal as the crystal-field splitting is progressively in- 
creased 1341. 
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FIG. 4. (Color online) Self-energies of the two-band model 
with increasing crystal-field splitting A. (a-b): Matsubara- 
frequency self-energy IniEm(J<^n) of (a) the lower band (fully 
filled for A > 0.5D), (b) the higher band (half-filled for 
A > 0.5D). Single-band (SB) results are shown with crosses, 
(c-d): Imaginary (c) and real (d) part of the real-frequency 
self-energy for the lower band (analytic continuation using 
stochastic maximum entropy). 



On Fig. |4l[a)(b) we display Iml]m(i(u;„) for imaginary 
(Matsubara) frequencies, for the lower and upper bands, 
respectively. The strength of correlations is characterized 
by the overall magnitude of I]„i(ia;„). As A increases, 
the occupancy of the upper band progressively dimin- 
ishes from ni = 1.5 (A = 0) to ni = 1 and the strength 
of correlations in the upper band increases, reflected also 
in Z^^ = 1 — d^{iu!)/d{iuj)\ui^Q corresponding to the 
enhancement of effective mass. As soon as the lower 
band is completely filled {712 = 2). which happens here at 
A = 0.5D, Si(za;„) saturates and becomes identical to 
the self-energy obtained in the single-band case A — s- 00 
(Fig.|4}D). This is in agreement with the qualitative argu- 
ments presented above: no particle-hole excitations can 
be created (at low energy) in the fully occupied band. 

Correspondingly, the correlations in the lower band 
progressively diminish, as the band is filled- in (Fig. |3Ji). 
However, in contrast to the upper band, no qualita- 
tive change is seen at A = 0.5D and the correlations 
evolve smoothly as A is further increased, remaining siz- 
able even at A = l.OD. On Figs. |4l^c)-(d) we display 
Iml]2(a; + J0+) and ReS]2(w + iO'^) analytically contin- 
ued on the real- frequency axis. At A = 0, the self energy 
has a pronounced peak at negative frequencies, corre- 
sponding to the pronounced lower Hubbard band. This 
structure persists through the regime A > 0.5 when the 
lower band is totally filled. To complete the picture, we 
display on Fig. [5] the energy-resolved spectral functions 
of the lower band. At the transition point A = 0.51?, the 
top of the lower band touches the Fermi level and devel- 
ops there a sharp quasi-particle peak at w = 0. Most of 
the spectral weight is carried, however, by a broad struc- 
ture seen at negative frequencies. When the lower band 
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FIG. 5. (Color online) Energy-resolved spec- 

tral function for the two-band model A2{e,ijj) = 
— Im [oj + ^ — e — A — E2(a; + iO^)] /tt. Left: contour 
map in the (e, to) plane, with —D < e < +D. Right: energy 
(aj)-resolved spectra for different energies e + A in the band. 
Upper panel: A = 0.5D. Lower panel: A = l.OD. 



lies below the Fermi level (A = l.OD), the quasiparticle 
structure merges with the broad Hubbard-like and the 
energy-resolved spectra are broad for most values of the 
momentum (except at the very top of the band). Hence, 
these model results fully confirm the qualitative expec- 
tations above and clearly demonstrate strong correlation 
effects in the fully occupied band due to inter-orbital in- 
teractions. 

In conclusion, we have shown that electronic correla- 
tions explain the spectral features seen in photoemission 
experiments on LaNiOa . Our DFT-I-DMFT results are in 
quantitative agreement with these measurements. While 
a description in terms of renormalized quasiparticles ap- 
plies at low-energy [35|, the 'kink' in the quasiparticle 



dispersion, shown here to be of electronic origin, signals 
the breakdown of this picture. We have demonstrated 
that inter-orbital interactions are responsible for the ob- 
served large broadening of the fully occupied t2g states. 
These correlation effects become important when the oc- 
cupied band is close to the Fermi energy and hence are 
relevant to other materials, for example to the broaden- 
ing of e' states in Na2;Co02 |36l438| . On a broader level, 
our results question the applicability of band theory to 
occupied bands strongly coupled to partially-filled ones. 
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